Theory for the nonlinear optical response of a nonspherical metal cluster 
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Using classical electrodynamics we determine the higher harmonic radiation by a nonspherical 
metal cluster in form of a matrix formalism. Extending the theory for the source of the higher 
harmonic radiation for spherical metal clusters as introduced by Ostling et al. [Z. Phys. D 28, 169 
(1993)] we calculate the sources for nonspherical particles. Employing these sources we develop the 
nonlinear Mie theory and determine the radiated higher harmonic fields generated by the cluster. 
Our theory is valid for arbitrary shape and arbitrary complex refractive index for cluster sizes much 
"j^Y , smaller and comparable to the wavelength of the incident light. 
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. The scattering behavior of clusters can be divided approximately into two regimes. First the scattering by clusters 
with sizes much smaller than the wavelength of light, the classical Rayleigh limit, which is dominated by microscopic 
effects and second the Mie range when the particle size becomes equal to or larger than the wavelength of light. The 
latter is characterized by resonances and interferences of surface electromagnetic modes described by classical physics. 
Since wc arc interested in this regime, we neglect microscopic phenomena. For spherical particles, the macroscopic 
effects known from linear Mie scattering |l]J| are even more pronounced in second and third harmonic generation, 
SHG and THG, respectively. Furthermore, they depend much more sensitively on the size and refractive index in the 
higher harmonic case For nonspherical particles, another important parameter governs the value and the angular 
dependence of the radiated intensities, viz. the curvature. Apart from the problem of determining this parameter as 
well as the size and refractive index from the angular dependence of the scattered intensities, the influence of clusters 
or well-defined nanostructures at surfaces on the SHG signal in reflection, especially the size dependence and the 
absolute value of the nonlinear response Q are of particular interest. An important advantage for the application 
of higher harmonic techniques is that they are, in contrast to linear optics, not influenced by background noise. 
Knowledge of the intensities radiated by nanoparticles deposited on surfaces or by surface roughnesses of nanometer 
size, should be one ingredient to determine the lateral spatial resolution of SHG theoretically, the lower limit of which 
is as low is as 1 nm || due to local field enhancements. Even in this range, the classical Rayleigh limit, an essential 
O ■ contribution results from the intensities radiated by nonspherical particles. 

O ■ In the literature one can find several approximate theories for the linear radiation from nonspherical clusters. The 
Rayleigh-Gans theory is valid for arbitrary shapes but only for materials where the absolute value of the complex 
refractive index is close to one (|e| ~ 1). Latimer |6| calculates the intensities scattered from ellipsoidal particles by 
relating every single point on the surface to an equivalent sphere and determines their linear response within Mie 
theory j7j. Other theories determine the intensities from aspherical particles of arbitrary size within a perturbation 
approach in terms of the linear deviations from the spherical shapes In the special case of spheroidal particles 
Asano and Yamamoto ||] use spheroidal coordinates which decouple the vector wave equations to obtain the exact 
results. 

Here we extend our theory for the nonlinear response of a spherical particle |]|ll|] to a particle of arbitrary shape. 
Therein the sources of the higher harmonic radiation are treated in terms of the surface charge a, which is, in the 
case of metals, equal to the normal projection of the polarization on the the surface (tr = n • P) and approximates the 
induced surface charge in the n-th harmonic case by the n-th power of the linear induced surface charge = (er^ 1 ))™ . 
This approach is motivated by the classical anharmonic oscillator model (see ) which is valid in the Mie regime and 
gives a good approximation for particles with can be considered as locally smooth. The n-th power of the surface 
charge, which is induced by the incident field acts as the source for the higher harmonic fields. It yields a discontinuity 
of the electrical displacement at the surface which can be expressed in terms of the boundary conditions for the radial 
component of the electrical fields. Using also the continuity of the tangential component of the electric field, the 
radiated fields can be fully determined. This is a simple task for spherical particles if spherical coordinates are used. 
In the case of nonspherical particles, however, matrix equations come into play where the size of the matrices depends 
on the ratio of the particle dimension to the wavelength of the light and also the geometry of the shape. 

To obtain the results for the linear scattering, which are necessary to determine the induced surface charge 
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a formalism by Barber and Yeh is used. This is based on an integral equation method by Waterman |L3 14 
and Schelkunoff's equivalence principle Jl5| yielding matrix equations. This is a good choice because Barber and 
Yeh use nearly the same representation of the fields with spherical coordinates and their theory is valid for arbitrary 
cluster-shapes and arbitrary complex indices of refraction. 



II. THEORY 

We determine the n-th harmonic electric field E^ radiated from arbitrarily shaped particles as a function of 
the n-th power of the linear surface charge cr^- 1 ' induced by the incident field which acts as the source of the n-th 
harmonic. The linear surface charge is given by = n • P' 1 ' where n is the vector normal at the surface of the 
cluster since in the case of metals only the normal component of the polarization is nonzero. The source a'™' = (a^) n 
induces a discontinuity of the normal component of the dielectrical displacement at the boundary of the particle. As 

(n) 

a consequence we use the conventional boundary conditions for the components of the electric fields to derive E^ u ( 
which then determines the scattering profile and all the characteristics of the nonlinear Mie scattering by nonspherical 
particles. Obviously 

n.(DW- D W)=W») 

and 

nxfEfi-E^O, 

together with D = eE, yield E^"| when is given. 

For the calculation all quantities, the fields and the source, are expanded in terms of spherical harmonics. This 
results in a matrix equation for the coefficient vectors of the fields and the source which completely determines 

(n) 

Eq U {.. We show in the following how this analysis can be performed in particular for the nonlinear Mie scattering of 
nonspherical particles, such an analysis had not been done before. 



A. Determination of the source a 



To determine the source we use the approximation a^ n ' = (cr*- 1 ^)". To obtain the n-th power of the linear 
charge we proceed as follows. The linear surface charge is defined as 

^(0 )¥ >) = n.pW(^) . (1) 

and can be expressed as 

a« (9, if) = i-Re [(1 - er 1 ) (E sc + E inc ) • n] e^* . (2) 

Here, we denote the scattered field by the index "sc" and the incident field by "inc" . Therein we refer to the boundary 
conditions of the normal component of the dielectrical displacement in the linear case n • (D sc + Di nc ) = n • Dj n (the 
index "in" states the field inside the cluster) and to the relation D = eE = E + 47rP. Furthermore, is expanded 
in terms of spherical harmonics 

v {1) {9, V ) = \ J2 4%(^) e " iUt+c - c - ( 3 ) 

l,m=±l 

From this it is easy to expand the n-th power of the linear surface charge in the form 

(9, V) = \Y, a tl Y ^ (0, <p) e" iTlwt + ex. . (4) 

l,m 

by raising the linear surface charge to the n-th power and neglecting terms with time-dependence different from 
e mut . Thus, we get the nonlinear sources <t'™'. As an example the coefficients of yield ( p,pT|) 
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= 2 E E a hAU / . (5) 
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i 1= l i 2 =l 

Therein the integrals over the Legendre Polynomials Yi iTn can be expressed in terms of the well known 3j-symbols. 
To calculate the coefficients the fields E out and E sc , which enter via Eq. |2|, are expanded in terms of vector spherical 
harmonics as introduced by Jackson jl0| in the form 

Ef } (x) = C (I) [K M (I, m) ft (hr) X,, m (9, <p) 

I, in 

777 [ 

+—K E {l,m)- r ^7xft(k l r)X l , m (e,ip) . (6) 
\m\ e(nu))k 



Therein, Xj im (0, <p) = LYj im (#, (p)/y/l(l + 1) is a vector spherical harmonic (L= l/i(r x V) is the angular momentum 

operator), C(l) = i l An(2l + 1), k — ntu/c and k\ — \J e(nu>)k. The multipole coefficients K^\l,m) and K^\l,m) 
refer to the magnetic (transverse electric TE) and electric (transverse magnetic TM) multipolcs. The index i specifies 
the fields external (i =out) or internal (i =in) to the cluster. The spherical Hankel functions f° nt (kr) — hi(kr) and 
Bessel functions f™(kr) — ji{k\r) describe the normal projection of the field inside and outside the particle. In the 
linear case the external field splits into the incident and the scattered field E inc and E sc , respectively. The coefficients 
of Ej, lc are known from the input field while E sc can be calculated by using the results from the theory for the linear 
problem by Barber and Yeh |Q . They expressed the coefficients of the internal field in terms of the coefficients of the 
incident field and then obtained the coefficients of the external field in terms of the coefficients of the internal field. 
In higher harmonic radiation where we will use the same expansion of the fields, however, only the external field 

(n) (n) 

Eq U ( radiated from the oscillating surface charge and the internal field E in occurs, since no incident field is present. 

B. Determination of the radiated field 

(n) 

To calculate the radiated field Eg Ut it is more convenient to introduce the abbreviations 
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K M (l,m) = C(l)K M {l,m) , K' E (l,m) = —^—C(l)K E (l : m) (7) 

| to I e{uj)k 

in the sum representation of Eq. |6[ This leads to 

E.[ ri) (x) = Y, [Km^^P (hr)X l>m (9,<p)+K' E (l,m)V x /,(*ir)X, im (0 j¥ >)] . (8) 

l,m 

In the n-th harmonic case the source acts as the discontinuity of the normal part of the electrical displacement 

tt .(DW-DW)=W*). (9) 

(according to our model for the nonlinear sources equals (a^ ) n ). Furthermore we will use the continuity condition 
of the tangential component of the electric fields 

n x ( E W - E^) = . (10) 

For particles of arbitrary shape the radius vector e r is no longer parallel to n and the radius r rather becomes a 
function of 9 and ip at the boundary since every point of the surface is determined by the angles. Furthermore the 
normal vector n now has the form 
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n(0, 4>) = n r (0, tp) ■ e r + n e (6>, • e g + n v (6>, <p) ■ e v . 

Here, (e r ,ee,e v ) are the basis vectors corresponding to spherical coordinates. As a result Eq. |l(] splits into three 
equations, one for every component 
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By using the series representations 
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V x fc, w (fcr)X, 



(12) 



multiplying every equation by Y/'. m ' where (V ,m!) runs over all (Z, m)-combinations used in the series representations 
of the fields and the source and integrating over the solid angle we get the following matrix equations 
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The coefficients are denoted by the vectors 
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where the coefficients in the vectors A^, and B^ are ordered in analogy to A^ 
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The matrix elements can 



be determined by using the expression for the curl terms in Eq. |12| with r depending on 9 and tp 
V x /,(fcr)X,, 



i a 

-3- [r/z(fcr)]e r x Xj . 
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g/Kfcr) /gr dY l;n 



^l(l + l)fi{kr)Y k , 

| 1 9r ayj, m 

sin 2 9V 3 9</3 



Then the matrixelements get the form 

Ml(l,m;l',m') = [ dtt (kr) 



BY,. 



1 dY, 



-n$- 
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dp 
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1 dr dY, 
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V^+T^^fcr)^, 
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sin 2 9 dip dip 
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Replacing the Hankel functions h\ x \kr) by the Bessel functions ji(kxr) the matrix elements for the matrices with the 



suffix 3 and 4 can be derived (see Eq. |12| ). Using as many (l',m') -combinations as (I, m)-combinations all matrices 
will be quadratic. Thus Eq. |l3J becomes a matrix equation off the form 
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which can be solved by standard numerical techniques. Moreover it is only necessary to calculate the coefficients A^. 

and A$ since they fully determine the n-th harmonic radiated field. This completes the solution of the nonlinear 
Mie scattering problem for arbitrary cluster shapes and arbitrary complex index of refraction. 

Note that the matrix equations and the many equations presented in this chapter are necessary to describe the 
interesting physics. They are a consequence of the nonsphericity of the particles implying boundary conditions that 
vary on the surface of the particle. This leads to the coupling of different electromagnetic modes. This is in contrast 
to the radiation by a spherical cluster. This feature is of general validity independent of the approximations made 
for the sources. One should study experimentally this conclusion of the theory. In the case of spherical particles all 
equations will collapse to those of the spherical case. 



III. DISCUSSION 

First we discuss the special case of a spherical cluster where no magnetic multipolcs are radiated and thus the 
coefficients A$ and are equal to zero. Furthermore the orthogonality of the spherical harmonics 

J Yi, m (9, 0) Yl m , (9, 4>) cM = S^S m . m , , (24) 
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can be applied in the integrals of Eq. |15|-|22| since the radius r and the normal vector n are independent of 9 and 
ip. Thus the matrices MJ and S 1,3 become equal to zero and the other matrices become diagonal. The coefficients 

A^(l, m) can then directly obtained from a["^. 

Thus deviations from the spherical shape cause that one mode of the source can excite several modes of the field 
which is expressed by the fact that the coefficients of the fields are coupled to the coefficients of the source by a 
system of equation. The shape of the particle directly governs the integrals in Eqs. |l5|-|22| and their nontrivial form 
for nonspherical particles leads to the system of equations in Eq. |i~3| which are valid for arbitrary cluster shapes. Also 
no restrictions to the size of the complex index of refraction are necessary in contrast to theories where a value of 
|e| « 1 is needed like the Rayleigh-Gans theory. 

Since our theory needs results from theories for the linear scattering it already includes all their numerical problems. 
In principle all modes must be summed in the series in Eq. ^ to represent nonspherical particles by spherical coordi- 
nates. So one limitation of the theory is the finite number of (I, m)-combinations which can be taken into account. 
As a consequence the matrices M and S in Eq. [l^ will become very large for strong deviations from the spherical 
shape. This purely numerical limitation is natural for particle shapes which cannot be described by coordinates which 
decouple the Helmholtz equation. 

On the other hand the theory has the advantage of a very compact form since the source for the nonlinear response 
is approximated just by <j(™) — (er^)" according to Ostling et al. [jll| and it only needs the boundary conditions for 
the electric fields from Eqs. || and [l(] for the derivation of the radiated fields from the source. The theory is especially 
valid for particles with sizes in the Mie range since it completely takes into account the combinations of the multipoles 
in the higher harmonic case (see Eq. ||) which are characteristic for the radiation of a particle in the Mie range. The 
approximation = (a^) n for the source term should be critically accessed if one is interested in details of the 
scattering profiles and effects due to the geometry of the cluster. 

One of the main purposes of this paper was to stimulate again work on Mie scattering from nonspherical particles. 
Clearly further work is necessary to demonstrate the validity of the analysis presented here. Our theory might 
stimulate the mathematical treatment of Mie scattering in the nonlinear case. 
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